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Abstract - Under a tilted magnetic field B, a two-dimensional electron gas (2DEG) with Rashba 
and Zeeman interactions is considered. The eigenvalues and their degeneracies, derived for various 
tilt angles, show crossing-free levels except at very high tilt. In this region, the tilt angle becomes 
the dominant factor in tuning the spin-splitting energy gap rather than the strength of the spin- 
orbit interaction. Moreover, concomitant with the crossings are the beats found in the 2DEG 
density of states. The crossings from different levels occur consecutively at around 87°. Similar 
new observations in Shubnikov-de Haas experimental measurements by Hatke et al. [5] attributed 
such phenomena to an in-plane magnetic field-induced increase in the effective mass. We show here 
that this behavior is inherent to a 2DEG where spin-orbit interaction and the in-plane magnetic 
field contribution are taken into account. 
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Introduction. — The advent of spintronics [T] and 
the quest to control spin in semiconductors have fueled the 
current research studies on Landau and spin energy lev- 
els of a two-dimensional electron gas (2DEG) PMj- Such 
investigations offer insight on how these energy levels and 
the spin-orbit coupling of electrons are affected by an ap- 
plied magnetic field and the 2DEG plane orientation. Tilt- 
ing of the 2DEG plane relative to the field direction was 
considered to enhance spin-orbit effects [5][7] making them 
easier to detect [5]. 

Early magnetotransport studies in Ga x Ini_ x As/InP [9 
at milliKelvin temperatures in strong magnetic fields at a 
tilt angle 9 of the order of ~ 80° showed coincidence of 
the first two Landau levels (LLs) with low carrier concen- 
tration, while the overlap disappears at sufficiently high 
mobility. This was believed to be caused by electron- 
electron interactions which shifts the system from a spin- 
unpolarized state to a spin-polarized state [5] . Other novel 
effects of tilting has since been observed which are not 
yet completely understood. For instance, the fractional 
quantum Hall states in the 2DEG in GaAs/AlGaAs het- 
erostructures was found to be suppressed as the sample's 
tilt angle 9 is increased, and resurging only at = 76° 
with a possible 2DEG nematic phase Moreover, ex- 
periments on the de Haas-van Alphen (dHvA) effect in 



high mobility 2DEG yield that magnetization oscillations 
with respect to an applied magnetic field exhibit beating 
patterns only at large tilt angles (9 > 75°) [JO]. The au- 
thors attribute these features to spin splitting induced by 
the structure inversion asymmetric potential in the het- 
erostructure. This type of coupling contribution is known 
as the Rashba spin-orbit interaction (SOI) [IT] . Usually 
the Rashba SOI is more dominant than the bulk-inversion- 
asymmetry in III-V [TJ] , II- VI [T3] and Si-based [TJ] semi- 
conductors. This offers a favorable opportunity for prac- 
tical applications because the Rashba spin-orbit coupling 
can be controlled via a gate voltage. 

In a different study carried out on a 2DEG near the sur- 
face of a highly doped p-type InSb crystal, measurements 
on the differential conductivity (dl/dV) also showed to 
have beating patterns of the LL intensity for different mag- 
netic fields [5] . Their results on the local density of states 
were consistent with the Rashba spin splitting. Not only is 
the Rashba SOI linked with the beats in these Shubnikov- 
de Haas (SdH) oscillations but it also gives way to a 
direction-dependent energy spectrum |15j . The authors 
of this present paper have shown previously that when 
both Rashba and Zeeman interactions are taken into ac- 
count there are no beating patterns in the density of states 
of the 2DEG for 9 < 80° gj. The authors observed that 
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Fig. 1: The eigenvalues as a function of the tilt angle for n 6 
[1,5]. A pair of identically colored lines indicates the two spin 
states of a given LL n. The dashed line in the inset points 
to the first crossing while that at the main frame the second. 
Here a = 3.5 x 10" 11 eVm and B = 0.1 T. 




a(10" U eVm) 



Fig. 2: The eigenvalues as a function of the Rashba SOI for 
n 6 [1, 5]. Here B — 0.1 T. A pair of identically colored lines 
indicate the two spin states of a given LL. The magenta lines 
serve as guides to the eye pointing to the crossing points. 



depending on the intensity of the applied magnetic field 
B, the Rashba SOI and Zeeman term signatures can be 
distinguished [5]. Moreover, at the tilt angles below 80°, 
the specific heat exhibited a suppressed beating pattern 
as a function of B. 

In this current work, we will distinguish the effects of 
9 from the rest of the tunable parameters. We will show 
how the energy spectrum of a two-dimensional electron 
gas changes as we sweep through the different directions 
of an applied magnetic field B. Special attention is de- 
voted to the large tilt case, that is when 9 90°, as it 
reveals the beats that seem to have only been restrained 
in smaller angles j4][5] . Another interesting feature consid- 
ered is that of a crossing which occurs when two adjacent 
LLs become degenerate for a given applied magnetic field, 
Rashba SOI strength, and tilt angle. We will show that 
these crossings are correlated to the beating patterns ob- 
served in thermodynamic quantities fTUlfTB] . 

Recently, SdH oscillations of ultrahigh mobility 2DEG 
in GaAs/AlGaAs quantum wells showed beating patterns 
at high tilt angles (> 87°) which indicate consecutive level 
crossings occurring at the same angle [5]. They explained 
their results as due to an increase in the carrier mass in- 
duced by the in-plane field with crossings dependent on 
the filling factor [6]. Here we will show that the density 
of states of a 2DEG model which incorporates Rashba 
SOI naturally exhibits multi-crossings at large tilt angles. 
From this a causal relationship of the LL crossings with 
the beating patterns observed in thermodynamic quanti- 
ties is established. 

Crossing States. — If Zeeman effect is the lone source 
of splitting, the SdH or dHvA effect is characterized by a 
single frequency [TTJIE] ■ This is because, just like the 
Zeeman term, the cyclotron energy huj c is also linearly 
proportional with the magnetic field magnitude B [17101?] , 



that is, uj c — eB±/m*. The perpendicular component B± 
is set along the z— axis and is equal to B z = B cos(6>) and 
m* is the effective mass. When the Rashba SOI is signif- 
icant, the spin-splitting is not simply linear with B. The 
Hamiltonian describing a single electron in this system can 
be written as 

h 2 P 

H = H a + H R + H z = h a(a X k) ■ z — ft - B, (1) 

2m* 

where Hq is the free particle energy, Hr is the Rashba 
spin-orbit interaction, and Hz is the Zeeman energy. 
Hence, on top of the strongly degenerate LLs, additional 
degeneracies can arise if the 2DEG further experiences a 
spin-orbit interaction. The strength of the Rashba SOI 
is indicated by the parameter a and is assumed to be 
constant in this work. Here a are the Pauli matrices 
and the magnitude of the wave vector k is determined 
by kj = — i\J j + j-Aj, where Aj is the j-th component of 
the vector potential. The Hamiltonian above is exactly 
solvable when B is normal to the 2DEG. However, if a 
finite in-plane component, B^ = Bsm(9), is present, the 
problem becomes undcrspecified and an analytic solution 
can only be achieved when additional specifications are 
imposed {41120]. 

We circumvent the underspecified problem by imposing 
equal probability among pairs of LLs that are known to 
cross when B\\ = [3J[S]. This yields an exact expression 
for the eigenvalues which is given as 

1 T 

E±=n + - + - (Vn + 1 + v/il) 
±\ v [ T (V^TT - v^) + 2ft,] 2 + 4ft+ft_. (2) 

Here n is the LL index, Qj = fJ-sBj /hu c is the j-th com- 
ponent of the rescaled Zeeman energy, T = 2a^/ (,/hui c is 
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the rescaled Rashba parameter, £ = m* /2H 2 and fl± = 
Q x ± iQ y . A detailed derivation of the eigenvalues can 
be found in 0]. In the results that follow, energy units 
are given in terms of the cyclotron energy fiuj c . In ad- 
dition, the 2DEG is considered to be on the x — y plane 
subject to an applied field in the x — z direction with mag- 
nitude B = \J B% + B\ . In general, m* is dependent on 
the energy due to nonparabolicity, but in this work we let 
m* — 0.05 x m e (m e is the free electron mass) which is 
the effective mass in a 2DEG in InGaAs heterostructures 
|21j and the same order in other semiconductor systems. 

In some theoretical studies H01122], the degeneracies at 
the crossing levels are lifted as 9 deviates from zero. That 
is not the case in this work for large tilt angles as illus- 
trated in fig. [TJ Although not included in the figure's 
scope, a crossing-free spectrum is observed in the region 
9 < 9 Xl where 9 Xl is the angle where the first crossing x\ 
shows up. As 9 is further increased more crossings begin 
to emerge. The crossing point x shifts to lower 9 as the LL 
index n is increased. This is expected since higher n cor- 
responds to a wider spin splitting as indicated in eq. ©■ 

Inspecting eq. we anticipate an increasing spin gap 
A with increasing 9. This gap between the two spin levels 
of a given LL — E~) expands until they cross their 
adjacent LL. The crossings, however, occur only at very 
large 9 as shown in the figure. 

It has been predicted that if a finite Bit is present, the 
energy levels do not cross as the Rashba SOI parameter 
is increased [3D]. We check the range of validity of this 
prediction regarding the energy levels as a function of the 
Rashba parameter, E(a), using our result in eq. ([2]). In- 
deed, for 9 < 6 Xl , E(a) is crossing- free [3]- However, for 
9 > 9 Xl , E{a) starts to exhibit crossings. This is depicted 
in fig. [21 At 88°, only one crossing is resolved because the 
gap widening is slow to change with respect to the increas- 
ing Rashba SOI strength. Relatively larger angles, 89° for 
example, exhibit a faster spread with a causing the oc- 
curence of multi-crossings. From the plot we understand 
that the tilt angle has a repulsive effect when 9 < 9 X \ but 
gives way to crossings when 9>9 X \. 

Comparing fig. \T\ and fig. &\ we can extract that large 9 
and small a both favor crossings. Increasing the Rashba 
parameter widens both the spin gap and the Landau gap. 
The interplay of these two produces the crossings at small 
a and the absence of crossings at large a as depicted in 
fig. [2J Whether Rashba SOI gives way or inhibits crossings 
is determined by 9. Note that in both figures, the Zeeman 
contribution is a constant for a given 9. Otherwise, it is 
diminished as the tilt increases. 

Figure [3] exhibits in detail how the spin gap varies with 
the Rashba SOI and the field's tilt angle. For the left 
panel, the A(0) curves relative to the tilt angle for two val- 
ues of a are qualitatively similar in shape. Only the mag- 
nitude is increased upon increasing a. However, at very 
large angles (9 > 85°) the two curves approach the same 
set of values. This signifies that, although the Rashba SOI 
strength brings about spin splitting of adjacent levels, the 
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Fig. 3: The spin splitting as a function of the tilt (left plot) 
and the Rashba parameter (right plot) when B = 0.1 T. 



large field tilting tunes the splitting behavior in addition 
to a. For the right panel, on the other hand, A(a) is linear 
with respect to a at small angles. At large angles, however, 
on top of its increased magnitude, the A(a) shifts from a 
linear behavior to a polynomial dependence on powers of 
a. Specifically at 9 = 88°, A has a cubic dependence on 
a. Figure 3 affirms how the primary control belongs to 9. 

Density of States. — In experiments, the density of 
states (DOS) are usually probed indirectly through mea- 
surements of transport properties such as the magnetiza- 
tion [TU] and the specific heat [23]. Recent developments 
have been able to probe the localized DOS in real space 
through scanning tunneling spectroscopy [2] [24] . The 
2DEG spatially averaged spectrum obtained in Ref. [2] has 
a beating pattern of the LL intensity in the first subband 
which the authors attributed to the Rashba spin splitting. 
Motivated by this observation, including the suppressed 
beating patterns in the specific heat at small and interme- 
diate angles [5] and the enhanced beats in the magnetiza- 
tion at large tilts [TUJ , the perceptible influence of 9 on the 
DOS in the presence of the Rashba SOI will be considered 
and shown here. This is done using the phenomenological 
form of DOS as defined by a series of Gaussian functions, 
namely, 



DOS(£) = 



eBj 



E 



1/2 



exp 



±\21 



2T 2 



(3) 

where 9 affects the eigenvalues E^ and the magnitude of 
B±. 

The broadening width T measures how the DOS devi- 
ates from its ideal impurity-free 2DEG behavior. Thus, 
this parameter is indicative of the number of impurities 
and/or the scale of disorder in the system. Extracted 
values of the broadening range from 0.4 - 1.5 meV in 
InAs/GaSb [H], 1-5 - 3.0 meV in InSb/InAlSb [25] and 
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Fig. 4: DOS as a function of the tilt angle at E ~ 16 meV 
for a = 3.5 x 1CP 1 eVm, B = 6 T, and different broadening 
widths. 

can reach as high as ^20 meV in InSb [3] depending on the 
acceptor concentration. The T values used here are of the 
same order as the theoretical fits in AlGaAs/GaAs which 
are dependent on the field magnitude, 0.1 meV x -JB(T). 
For simplicity, we keep T constant. The case of T < 1.0 
meV ensures that the spin-splitting and the tilt effects are 
not washed out and dominated by the broadening. 

As shown in the inset of fig. the DOS peaks 
monotonously decreased for 9 < 86°. This behavior can 
be traced to the direct proportionality of the prefactor of 
eq. © to B±. When 9 is large as in the case of 9 > 86°, 
the vanishing amplitude at the beat node suggests that the 
two spin subbands oscillate with similar frequencies. As 
9 is decreased, the difference between the frequencies in- 
creases. When 9 < 86°, for example, the two frequencies 
become very different that only one of them dominates. 
This gives rise to the single- frequency-DOS. In Rcf. [15] , 
the single-frequency oscillation was observed in the high 
B region of the spectrum while the two-frequency beating 
pattern was detected in the low B region. This is anal- 
ogous to our observation where the latter corresponds to 
the high 9 domain. Furthermore, we observe that the DOS 
oscillations with double frequencies, namely the beats, be- 
come more pronounced as B is intensified. This is due to 
the proportional increase in DOS amplitude. 

Referring still to fig. SI we now relate the beats with 
the crossings. Both beats and crossings are absent until 
9 reaches a certain large value (> 86°). We deduce that 
each DOS antinode is concomitant with a crossing point. 
Since two energy levels merge to the same energy value E, 
the degeneracy of that E is doubled. This translates to 
a maximum DOS envelope function. The crossings occur 
when A = jhuj c , where j = 1, 2, 3, ... On the other hand, 
nodes result when the LLs are equidistant to their adjacent 
levels. Such a condition is satisfied when A = (j/2)Tiio c , 
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Fig. 5: DOS versus B± at E ss 16 meV showing the phase shift, 
the spin splitting and the crossings. The curves are shifted for 
clarity. Here a = 3.5 x 1CT 11 eVm and T = 0.5 meV. The red 
arrows point to the spin splitting of a LL while the blue arrows 
to the reemergence of the spin levels after crossing. 

where j here corresponds to the order of the node [T7] , 
The latter resembles the usual evenly spaced harmonic 
oscillator spectrum. 

In fig. 2J we can observe that a large T can mask the 
presence of the beats. This might help explain the absence 
of SdH beats in materials that are known to have strong 
Rashba SOI [55] ■ In addition, when the tilt angle has not 
yet reached the threshold value 9 X \ where the beats begin 
to emerge, a finite B\\ only suppresses the beats and the 
SdH is seen to oscillate with a monotonously decreasing 
amplitude. In this case, for all 9 < 9 X \ the beats are 
absent. 

Collating all of the abovementioned observations, beats 
are found at the two ends of the tilt spectrum - the zero 
tilt {B\\ = 0) [J and the large tilt (B\\ -> B). In be- 
tween these two extremities, the beats are not observed. 
Moreover, only at these regions did we also observe the 
crossings which is indicative of beats and crossings being 
concomitant. Since the beats are due to the spin-splitting 
phenomena induced by the Rashba SOI and crossings are 
found in the case of no tilt [2J[35] but are not found for 
9 < 80° [1], then the magnetic field's tilt relative to the 
2DEG plane shifts the spin splitting-induced crossings at 
large tilt as observed in this paper. 

To quantify how large the large-tilt region is, we plot- 
ted DOS verus B± = B z . See fig. [5] We found out that 
the DOS versus B± has characteristic 180° phase rever- 
sals at 9 rl w 84.75°, 9 r2 w 88.5° and 9 r3 w 89°. The DOS 
are in-phase for < 9 < 6 r i, for 9 r \ < 9 < 9 r 2, and for 
9 r 2 < 9 < 9 r 3. The complete phase reversal, with the same 
frequency, however, is limited only to the small B± region 
(B± < 1 T) where the LL splitting is negligible. This has 
been observed experimentally [2JJ . When the splitting be- 
comes evident, it is not straightforward to compare phases 
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because the LLs begin to split into two spin levels (See red 
arrows.). As 9 is increased, the gap between the spin lev- 
els widens until they cross other LLs (See thickened line 
at 9 = 87°.). Compared to 6 = 86°, each amplitude at 
9 = 87° is approximately doubled and the widths of the 
peaks became narrower. These arc evidences of the cross- 
ing of two spin levels of neighboring LLs. Note that here 
we observe multi-level crossings at 9 = 87°, specifically 
at B± = 1.126 T, 1.408 T and 1.879 T. This is in agree- 
ment with the experimental data of Hatkc et ai, in which 
this unusual observation of level crossings occurring at the 
same angle are attributed to an in-plane magnetic field- 
induced increase of the carrier mass. Indeed in Ref. [35] 
the distortion on a Fermi loop was directly measured us- 
ing the magnetic electron focusing effect in a 2DEG. This, 
they say, affects the effective mass because of the antisym- 
metry of the cyclotron orbit. Although, we kept m* in our 
work a constant for the simulated data presented herein, 
the eigenvalues in eq. ^ are derived from first principles. 
The exact derivation incorporates the Rashba SOI and the 
in-plane field whose mutual influence on the electron's or- 
bit changes the characteristics of the energy spectrum and 
consequently, the density of states. 

It is important to mention that we do not observe here 
the anticrossings of adjacent spin-split LLs seen in Ref. [7]. 
They attribute the anticrossing to the spin-orbit coupling 
of the collective spin configuration of the quantum hall fer- 
romagnets and to the nonparabolicity causing spin mixing. 
The authors of [53] believe that it does not come from 
electron-electron interactions but from the pronounced 
nonparobolicity of the bands. From our work, we infer 
that anticrossings must come from nonparabolicity of the 
energy bands that arc induced by other interactions which 
are not from the structure inversion asymmetric potential 
and the tilting of the applied field. 

In conclusion, we studied how Landau level crossings 
emerge at large tilting. These crossings bring about beats 
in the density of states which are manifestations of two 
spin subbands oscillating at similar frequencies. In agree- 
ment with recent experimental data from Hatke et ai, 
the crossings are found to occur consecutively at around 
87°. While they attribute this unusual behavior to the bi- 
plane magnetic field-induced carrier mass enhancement, 
this work has shown that the multiple crossings are in- 
trinsic to a two-dimensional electron gas with Rashba spin- 
orbit interaction and under intense in-plane magnetic field. 
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